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ABSTRACT 


The microlensing surveys, such as OGLE or MACHO, have led to the discovery of thou¬ 
sands of RR Lyrae stars in the Galactic Bulge and in the Magellanic Clouds, allowing for de¬ 
tailed investigation of these stars, specially the still mysterious Blazhko phenomenon. Higher 
incidence rate of Blazhko (BL) variables in the more metal-rich Galactic Bulge than in the 
LMC, suggests that occurance of Blazhko effect correlates with metallicity (Moskalik and 
Poretti 2003). To investigate this problem, we calibrate the photometric method of determin¬ 
ing the metallicity of the RRab star (Kovacs and Zsoldos 1995) to the /-band and apply it to 
the OGLE Galactic Bulge and LMC data. In both systems, metallicities of non Blazhko and 
Blazhko variables are close to each other. The LMC Blazhko pulsators prefer slightly lower 
metallicities. The different metallicities of the Galactic Bulge and the LMC, can’t explain the 
observed incidence rates. 

As a by-product of our metallicity estimates, we investigate the luminosity-metallicity 
relation, finding a steep dependence of the luminosity on [Fe/H]. 

stars: horizontal branch - stars: oscillations - stars: variables: RR Lyr - stars: abundances 


1 Introduction 

The Blazhko effect, discovered for the first time in 1907 for an RR Lyrae star 
RW Draconis (Blazhko 1907), appears as a cyclic modulation of shape and 
amplitude of the light curve. About 30% of all galactic, fundamental mode 
RR Lyraes (RRab stars) display this effect (Smith 1995). Despite almost a 
century of study, we still don’t understand this phenomenon. 

The microlensing surveys, such as the Optical Gravitational Tensing Exper¬ 
iment (OGLE; Udalski et al. 1992) or MAGHO (Alcock et al. 1996), have led 
recently to the discovery of thousands of RRab stars in the Galactic Bulge and 
in the Large Magellanic Gloud (LMC), hundreds of them displaying the Blazhko 
effect. Detailed analysis of the MACHO LMC data (Alcock et al. 2000, 2003) 
showed, that Blazhko variables are multiperiodic pulsators with several closely 
spaced frequencies. The primary peak of the spectrum, /o, is accompanied 
by close secondary frequency (dublet) or by two secondary frequencies located 
symmetrically on both sides of the primary peak (symmetric triplet). This 
structure is repeated at the harmonics of the primary frequency. There are also 
nonequidistant triplets and multiplets present in the spectra, but these are rare 
cases (Alcock et al. 2003). The separation between frequencies. A/, is usually 
very small, typically A/ < O.lc/d. Corresponding Blazhko periods, Pb = 1/A/, 
are of order of days to hundreds of days. Eollowing Alcock et al. (2003), we 
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Table 1: Statistics of RRab stars in the Galactic Bulge and in the LMC. Ns 
and Nbl denote the number of RRO-S and RRO-BL stars respectively. 


Object 

Catalogue 

Ns 

Nbl 

% 

Ref. 

Galactic 

OGLE-I 

112 

35 

23.4% 

Moskalik, P., Poretti, E. 2003 

Bulge 

OGLE-II 

1556 

472 

23.3% 

Mizerski, T. 2003 

LMC 

OGLE-II 

4608 

843 

15.5% 

Soszyhski, 1. et al. 2003 


MACHO 

4882 

751 

13.3% 

Alcock, C. et al. 2003 


denote singly periodic, fundamental mode variables as RRO-S, variables with 
dublets as RRO-BLl and variables with symmetric triplets as RR0-BL2. The 
earlier analysis (Alcock et al. 2003) showed, that latter two types overlap, in¬ 
dicating the same underlying phenomenon, so the common notation, RRO-BL, 
will also be used. 

Statistics on recently analyzed objects is given in Table 1. 

It is clearly visible that Blazhko phenomenon is more frequent in the Galactic 
Bulge, compared to LMC. Moskalik and Poretti (2003) suggested, that this 
is due to the difference in metallicity, which is higher in the Galactic Bulge. 
Verification of this hypothesis, which is the goal of this work, requires an effective 
method of estimating the metallicity of RRab star. Such a method was proposed 
by Kovacs and Zsoldos (1995) for RRO-S variables. Metallicity is determinated 
on the basis of the Fourier decomposition parameters of the F-band light curve. 
The OGLE data, which is the subject of this analysis, is mostly the Cousins 
/-band photometry, so the new calibration of this method is necessary. 

Information about the Blazhko pulsators metallicity may be also interest¬ 
ing for theoretical modeling. Currently considered models incorporate nonra- 
dial geometry of pulsation {e.g., Nowakowski and Dziembowski 2001, Shibahashi 
2000), but none can reproduce all of the observed features. Any information 
about [Fe/H] may be very useful. 

In Section 2 we describe in detail the method of estimating iron abundance 
of the RRab star, its calibration to the /-band and application to Blazhko vari¬ 
ables. In Section 3 the techniques of light curve analysis both for singly peri¬ 
odic and Blazhko variables are described. In Section 4 we analyze the Fourier 
parameters of the light curves and metallicities of the RRab variables from the 
OGLE-I Galactic Bulge sample and the OGLE-II LMC sample. The luminosity- 
metallicity relation for the LMC RRO-S pulsators is examined. Conclusions are 
collected in Section 5. 
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2 Method 

2.1 The photometric metallicities 

Following the idea, that shape of the light curve depends on physicall parameters 
of the star, Kovacs and Zsoldos (1995) proposed a new method of estimating 
metallicity of the RRab star, based on its light curve Fourier parameters. The 
method was further developed by Jurcsik and Kovacs (1996, hereafter JK). They 
found, that the basic relation between [Fe/H] and Fourier parameters is linear 
and contains only the period and the Fourier phase, Their method was 

calibrated on 84 singly periodic, galactic field RRab stars, with good V light 
curves and spectroscopic metallicities. All the [Fe/H] values were transformed 
to a common metallicity scale introduced by Jurcsik (1995) and based on the 
recent High Dispersion Spectroscopy (HDS) measurments (hereafter refered as 
the HDS scale). This scale provides consistent metallicities for both field and 
cluster RR Lyrae stars. Jurcsik (1995) also gives a new, general AS [Fe/H] 
transformation formula. 

The method was tested by JK for 11 Galactic and LMC globular clusters. 
Photometric metallicities, averaged for several RR Lyraes, ([Fe/Hjgt), were com¬ 
pared to cluster’s spectroscopic metallicity, [Fe/H]obs- It was found, that photo¬ 
metric abundances generally fit well spectroscopic values, but in the low metal¬ 
licity range, [Fe/H]obs<— 2.0, photometric abundances are higher, on avarege by 
(([Fe/H]fit) — [Fe/H]obs) = 0.34±0.04 (3 clusters). No explanation was given for 
this discrepancy. In higher metallicity range, [Fe/H]obs > —1.8, the discrepancy 
still exsist, but it is very small, namely (([Fe/Hjgt) — [Fe/H]obs) =0.05 ±0.03 (8 
clusters). Since the HDS scale is not strictly reproduced by photometric method, 
we will refer to the photometric metallicities obtained from JK formula as based 
on the JK photometric scale. 

2.2 Calibration in the /-band 

As a by-product of the OGLE project, a good quality photometry for thou¬ 
sands of RRab variables has been collected. Since this photometry is mainly 
in the standard Gousins /-band, estimating metallicity of the star requires new 
calibration of the photometric method. Unfortunately there is much less good 
photometry for field RR Lyrae stars in the /-band, than in the V, so the cali¬ 
bration is based only on 28 singly periodic, fundamental mode pulsators*. The 
HDS metallicity scale is adopted here. 

For 21 field RRab stars good / light curves are available. For three of them, 
DX Del, V445 Oph and VY Ser, photometry is in the Johnson’s /-band. Trans¬ 
formation to the standard Gousins / is done, using Bessel (1983) formulae. Ex¬ 
cept V Ind, all the stars are also in JK calibrating sample, so their metallicities 

•Comparing the Clementini et al. (1990) V photometry, with Lub (1977), originaly Wal- 
raven photometry, JK claimed, that V Ind shows the Blazhko effect. However extensive 
Clementini et al. (1990) photometry alone, shows no sign of Blazhko effect, so following these 
authors, we treat V Ind as a single mode pulsator. 



4 


Table 2: The calibrating data set. Globular cluster star Id’s are consistent with 
Sawyer-Hogg (Hogg 1973). Light curve references: (1) Liu, T., Janes, K.A. 
1989; (2) Skillen, 1. et al. 1993; (3) Hansen, L., Petersen, J.O. 1991; (4) Barnes, 
T.G. et al. 1988; (5) Glementini, G. et al. 1990; (6) Fernley, J.A. et al. 1990; (7) 
Cacciari, C. et al. 1987; (8) Mazur, B. et al. 1999; (9) Kaluzny, J. et al. 1997. 


Star 

Ref. 

[Fe/H] 

Star 

Ref. 

[Fe/H] 

SW And 

1 

-0.06 

AV Peg 

1 

-f0.08 

WY Ant 

2 

-1.39 

AR Per 

1 

-0.14 

V Cae 

3 

-1.71 

BB Pup 

2 

-0.35 

RR Get 

1 

-1.29 

VY Ser 

6 

-1.58 

W Crt 

2 

-0.45 

V440 Sgr 

7 

-1.21 

DX Del 

4 

-0.32 

W Tuc 

5 

-1.37 

SU Dra 

1 

-1.56 

TU UMa 

1 

-1.15 

RX Eri 

1 

-1.07 

NGC 6362 vl9 

8 

-0.83 

V Ind 

5 

-1.24 

to Gen v40 

9 

-1.33 

RR Leo 

1 

-1.30 

iv Gen v44 

9 

-1.14 

SS Leo 

6 

-1.56 

to Gen v51 

9 

-1.37 

TT Lyn 

1 

-1.50 

iv Gen v62 

9 

-1.35 

RV Oct 

2 

-1.08 

to Gen vl02 

9 

-1.56 

V445 Oph 

4 

-bO.Ol 

CO Gen vll5 

9 

-1.59 


are adopted. Metallicity for V Ind was obtained by Layden (1994) and transfor¬ 
mation to HDS scale is done according to formulae given in JK (their equation 
2 ). 

Photometry for another six variables comes from u Gentauri (Kaluzny et al. 
1997), which is the globular cluster with high metallicity spread. Photometry 
is given in the instrumental system, but according to Kaluzny et al. (1996) it 
reproduces the Cousins /-band fairly well, and transformation is not necessary. 
Metallicities of these variables were obtained using Caby photometry technique 
(Rey et al. 2000), calibrated on Layden (1994) stars, so the transformation to 
the HDS scale is done just as above, using JK transformation. 

The last variable comes from NGC 6362 (Mazur et al. 1999), with metallicity 
obtained by Costar and Smith (1988) and transformed to the HDS scale by 
Jurcsik (1995). 

All objects with light curve reference and [Fe/H] value are presented in 
Table 2. 

Each light curve is fitted with the Fourier series of the form: 

N 

I{t)=Ao + '^AkSin{2TTkft + (j)k) ( 1 ) 

fe=i 


The order of the decomposition, N, is chosen to satisfy the condition: AkfaAk > 
4, for all the amplitudes. aA^ stands for the standard error of the amplitude, 
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Ak- The pulsation frequency, /, is adjusted to obtain the minimal value of the 
fit’s dispersion, a. In the case of WY Ant, optimal period is strongly dependent 
on the decomposition order, so the Kholopov (1995) period is adopted. Periods 
of w Centauri variables are derived from the more extensive V photometry. The 
TU UMa period is obtained by phasing the data from two sources (Liu and 
Janes 1989, Barnes et al. 1988). In all cases, points deviating from the fit by 
more than 4(T are rejected. 

Linear combinations of the Fourier parameters, up to the fifth order, namely 
the period, P, the amplitudes, Ai—A^, and the Fourier phases, (j)2i—4>5i (</>fci = 
(j)k — k<j)i), are tested to fit the observed metallicities, [Fe/H]obs- The best two- 
parameter formula contains, like in JK, the period and the <j)^i phase: 

[Fe/H]= -3.142 -4.902P 

±0.646 ±0.375 ±0.104 cr = 0.18 

a denotes dispersion of the fit. Abundances determined from above equation will 
be denoted by [Fe/H] 2 . The three-parameter formulae have smaller dispersions. 
The best one, contains additionaly amplitude A 2 : 

[Fe/H]= -6.125 -4.795P ±1.181(()3i ±7.876^2 , . 

±0.832 ±0.285 ±0.113 ±1.776 cr = 0.14 

Abundances determined from above equation will be denoted by [Fe/Hjg. For¬ 
mulae using four or more parameters are much worse. 

Observed vs. fitted abundances for the calibrating sample are presented on 
Fig. 1. 

2.3 Application to Blazhko variables 

Light curve shape of the singly periodic RRab variables, used in the calibration, 
is time-independent. In the case of the Blazhko variables, shape of the light curve 
varies from season to season. Light curve phased with the fundamental period 
is characterized by a large scatter, specially in the light maximum (Fig. 2a). 
Decomposition of such a curve is useless. Jurcsik, Benkd and Szeidl (2002) 
showed, that also seasonal curves, at different phases of the Blazhko cycle, don’t 
fit the typical singly periodic curve. However, we can extract the fundamental 
mode light curve, by prewithening the data with the additional close frequencies 
(Fig. 2b). We assume, that Fourier parameters of this light curve may be used 
to estimate the Blazhko variable metallicity. 

To prove this assumption, one should have some Blazhko variables with spec¬ 
troscopic metallicities, and extensive photometry, allowing for detailed Fourier 
analysis and filtering out additional close components, in order to estimate pho¬ 
tometric metallicity for comparison. Such a /-band photometry is available for 
the Galactic Bulge and the Magellanic Clouds OGLE samples, but there is only 
one individual metallicity value for Blazhko variable of the LMC, that is also in 
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our sample, namely SC_21 16249 (OGLE Id)^. Spectroscopic metallicity value 
comes from Gratton et al. (2004), and is transformed to HDS scale as described 
in Section 4.2, leading to value —1.63±0.11. Estimated, photometric metallicity 
(using Eq. 2, cf. Sec. 4.2) is —1.35±0.18, so the difference between these values 
is —0.28±0.21 - quite large, but still acceptable (1.3 ct). 

Of course one star is unsufRcient to test the method, however, for some galac¬ 
tic field, Blazhko pulsators, extensive and homogeneous V photometry suitable 
for Fourier analysis exists, and can be used to test the method in case of the 
V band. Unfortunately some of them (V442 Her, Schmidt and Lee 2000; GV 
And, Lee et al. 2002) don’t have spectroscopic [Fe/H], or value is very uncertain 
(AH Cam, Smith et al. 1994). In the case of XZ Dra, and RS Boo, available 
homogeneus sets of V photoelectric measurments (Szeidl et al. 2001, Kanyo 
1986) were made for 0-C analysis, so only the rising branch and maximum of 
light are well covered. We are left with only five stars, for which we are able to 
test the method, namely: DR And (Lee and Schmidt 2001), XZ Cyg (LaCluyze 
et al. 2004), AR Her (Smith et al. 1999), RR Lyr (Smith et al. 2003) and RV 
UMa (Kovacs 1995). For RV UMa detailed Fourier analysis of Flight curve was 
already done by Kovacs (1995), so we use their extracted Fourier parameters 
of the fundamental mode light curve (their Table 3) to calculate the photomet¬ 
ric metallicity. In case of DR And, XZ Cyg, AR Her and RR Lyr we perform 
the same analysis as described in Section 3 (and in greater detail in Appendix) 
to extract Fourier paramameters of the fundamental mode light curve. Calcu¬ 
lated photometric metallicities (JK’s Eq. 3) are presented in column 3 of Table 
3. Their errors are assumed to be 0.15 dex (JK, Jurcsik 1998 and references 
therein). 

In case of RR Lyr, XZ Cyg and RV UMa, spectroscopic metallicities are 
weighted averages of the Lay den (1994) and Suntzeff (1994) values, transformed 
to the HDS scale as described in JK. We used the same weightening scheme 
as JK, to set the metallicities on the same scale, as the V photometric method 
was calibrated. Metallicity for DR And comes from Layden (1994) and is trans¬ 
formed as in JK. In case of AR Her spectroscopic [Fe/H] is derived from Preston 
(1959) AS, using transformation given by Jurcsik (1995). Spectroscopic [Fe/H] 
values are collected in column 2 of Table 3 together with errors of the weighted 
averages. In case of AR Her the error coresponds to the uncertainity of one unit 
in AS. Difference between photometric and spectroscopic values is presented in 
column 4. 

In case of DR And, XZ Cyg, AR Her and RR Lyr there is a nice agreement 
between estimated and observed [Fe/H]. Only in case of RV UMa the difference 
is larger, but acceptable (1.2 ct). 

We conclude that the described method may be used to obtain reliable values 
of Blazhko variable metallicity in case of the V light curve and assume that it 
also holds for the I band, which is used through this paper. 

tWe mention that also Borissova et al. (2004) gives individual metallicity values for some 
LMC RRab stars, six of them being Blazhko variables (Alcock et al. 2003). However they 
mix the spectroscopic values with photometric metallicities (obtained from V light curves) 
and hence we can’t use their values, since we do not know, which ones are spectroscopic. 
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Table 3: Spectroscopic vs. F-photometric metallicities for Blazhko variables. 


Star 

[Fe/H]spec 

[Fe/H]phot 

A[Fe/H] 

DR And 

-I.22±0.07 

-l.I4±0.I5 

-0.08±0.17 

XZ Cyg 

-I.18±0.05 

-1.07±0.I5 

-0.11±0.16 

AR Her 

-I.ll±0.19 

-l.I4±0.I5 

0.03 ±0.24 

RR Lyr 

-I.12±0.01 

-l.I0±0.I5 

-0.02T0.15 

RV UMa 

-0.90±0.11 

-l.I3±0.I5 

0.23T0.I9 


3 The data and techniques of light curve analy¬ 
sis 

Singly periodic variables. The OGLE-I catalogue (Udalski et al. 1994, 1995a, 
1995b, 1996, 1997), contains 112 RRO-S stars. The OGLE-II EMC sample 
(Soszyhski et al. 2003), is almost 40 times more numerous, and contains 4607 
singly periodic, fundamental mode pulsators. Procedure of the analysis is the 
same for both samples. Light curves are fitted with the Fourier series, with 
period adjustment and bad data points rejection, just as described in Section 2.2. 
Since the data covers multiple seasons, possible instrumental drifts are modeled 
by a sine term with period of 50,000 days added to Eq. (1). 

From the Galactic Bulge sample we exclude MM7_V63, which is probably 
a member of the Sagittarius Dwarf Galaxy (Olech 1997). In case of the LMC 
sample, about 200 stars have noticably higher apparent luminosities, than the 
others. Those which are half a magnitude brighter than the sample average for 
a given period, are treated as blends and excluded from the sample. Only a 
good quality curves, that are decomposed to the third or higher order are used 
in later analysis. This is required by the method of estimating metallicity of 
the star. 90 RRO-S variables from the Galactic Bulge and 3990 RRO-S from the 
LMC passed the described criterions. 

Blazhko variables. In order to estimate the Blazhko variable abundance, 
the Fourier parameters of the fundamental mode must be extracted. This is 
done by fitting the light curve with the Fourier series, containing all the fre¬ 
quencies, also that corresponding to the additional close components. Then, 
the parameters corresponding to the primary frequency and its harmonics, are 
chosen to the further analysis. 

The Galactic Bulge sample was already studied in detail by Moskalik and 
Poretti (2003). They gave accurate frequencies and all the Fourier terms present 
in the spectra for 35 Blazhko variables present in the OGLE-I sample. We just 
use their Fourier series form and choose the proper parameters to the further 
analysis. Five stars with more complicated features (not being simple BLl or 
BL2 type, but for example, nonequidistant triplet, variable frequency values), 
are also analyzed. 

The case of the LMC data is more complicated. Soszyhski et al. (2003) pub- 



lished only the list of multiperiodic variables, with primary frequency, and one 
additional close frequency of the highest amplitude, not distinguishing between 
RRO-BLl and RR0-BL2 types. Extraction of the primary mode parameters 
is done by an automatic procedure, that searches only for BLl or BL2 types. 
Details of the procedure are given in the Appendix together with a test of the 
algorithm based on known results from the Galactic Bulge Blazhko pulsators. 

Not counting the blends (excluded on the basis of the same criterion as 
above), 702 LMC stars passed the algorithm criterions, 58% being of type BLl, 
and 42% of type BL2. Corresponding numbers in the MACHO sample (Alcock 
et al. 2003) are 55% (BLl) and 45% (BL2), so the results are consistent. 

4 Results 

4.1 The Galactic Bulge pulsators 

The singly periodic RR Lyrae pulsators from the OGLE-I sample, were already 
analyzed by Poretti (2001). Moskalik and Poretti (2003) considered the Blazhko 
pulsators. 

Eor all the stars, abundances are estimated using both two and three-parameter 
formulae. Metallicity distributions are presented on Pig. 3. In the case of singly 
periodic pulsators both formulae (Eq. 2, 3) give consistent results. Averaged 
metallicities are ([Fe/H] 2 ) = —1.04±0.03 and ([Fe/Hja) =—1.05±0.03. 

Estimated abundances are contained approximately in a range (—1.4,—0.8) 
for the RRO-S stars. Since the accuracy of the two-parameter formula is only 
(7 = 0.18, it is necessary to check, whether the spread of estimated [Fe/H] is 
physically significant. We do this by constructing the period-amplitude diagram 
(the Bailey diagram). It is a well known fact (Smith 1995), that for a given 
period, the RRab stars of lower metallicity have higher amplitudes. As we can 
see on Fig. 4, using metallicities obtained with our formulae, we reproduce this 
relation satisfactorily. Dispersions of metallicities (Fig. 3), larger than the error 
of the method, indicate intrinsic metallicity spread, of the order of ~0.16 (using 
two-parameter formula). 

Both of our [Fe/H] estimating formulae contains Fourier phase, (/ 31 . It is 
interesting to see where do stars of different [Fe/H], fall in the period-phase 
diagram. Plots of (psi and (f> 4 i vs. period for RRO-S stars, are presented on 
Fig. 5. Of course in case of two-parameter formula (left panel of Fig. 5) and (/31 
plot, lines of equal [Fe/H] are straight. The main sequence is formed by stars of 
[Fe/H] close to the sample average. The area above and below this sequence, is 
occupied by stars of higher and lower [Fe/H] respectively. This picture remains 
true for higher order phases as can be seen on (/41 plot. Also addition of A 2 
amplitude to the [Fe/H] estimating formula (right panel of Fig. 5), doesn’t 
change this picture significantly. Connection between location of the star on the 
period-phase diagram and its metallicity, for the same Galactic Bulge sample, 
was already noticed by Poretti (2001). The author distinguished three discrete 
tails, which he connected with different [Fe/H]. On the base of our study, we 



9 


Table 4: Averaged Fourier amplitudes and phases in three period bins, both for 
singly periodic and Blazhko variables from the Galactic Bulge. The error of the 
mean is given in a parenthesis. 



P<0.5 

Pe (0.5,0.6) 

P>0.6 


RRO-S 

RRO-BL 

RRO-S 

RRO-BL 

RRO-S 

RRO-BL 

Ai 

0.251(5) 

0.227(5) 

0.203(6) 

0.20(1) 

0.162(9) 

0.13(2) 

A2 

0.128(3) 

0.097(6) 

0.106(4) 

0.093(7) 

0.076(6) 

0.06(1) 

A3 

0.091(3) 

0.051(6) 

0.070(3) 

0.057(5) 

0.045(5) 

0.029(6) 

A4 

0.059(2) 

0.035(4) 

0.044(2) 

0.033(4) 

0.026(3) 

0.021(4) 

As 

0.040(2) 

0.024(4) 

0.028(2) 

0.025(2) 

0.018(2) 

0.013(-) 

4>21 

2.57(2) 

2.59(3) 

2.78(3) 

2.84(5) 

3.12(3) 

3.14(8) 

031 

5.37(5) 

5.25(9) 

5.83(4) 

5.8(1) 

6.39(7) 

6.4(1) 

041 

1.92(7) 

2.3(6) 

2.58(7) 

2.9(4) 

3.4(1) 

3.4(3) 

051 

4.73(9) 

4.3(4) 

5.59(9) 

5.0(1) 

6.5(2) 

6.7(-) 


may state, that better description incorporates areas of different [Fe/H], that 
are parallel to the sequence formed by most of the system stars, with continuous 
transition between them. 

Now, let us compare some features of the non-Blazhko and Blazhko variables 
light curves. On Figs. 6 and 7, the Fourier amplitude of the fundamental mode, 
Ai, and higher order amplitudes, A 2 — A 5 , and phases, (/> 2 i—</' 5 i, are plotted 
against period. The pulsation periods range from 0.3 to 0.9 day, being shorter 
than 0.7 days for Blazhko pulsators. On the plots of Fourier phases and Ai 
vs. period, there is no clear distinction between RRO-S and RRO-BL variables. 
However on the plots of A 2 — A 5 , it is visible, that BL variables prefer lower 
values of amplitudes. The effect is clearly visible, if average values, for different 
period bins, are calculated. Results are presented in Table 4. For all period 
bins, all amplitudes are smaller for BL pulsators. The effect is weakest for 
Ai. The plot of amplitudes A 2 — A^, vs. the amplitude, Ai , is presented on 
Fig. 8 . For low Ai values, Blazhko pulsators and singly periodic stars follow the 
same sequence, but for Ai > 0.2, Blazhko pulsators prefer lower amplitudes of 
harmonics. 

Noticed relations among Fourier amplitudes are reflected in metallicity dis¬ 
tributions for Blazhko pulsators (Fig. 3). The two-parameter formula gives the 
similar shape of the distribution as for the singly periodic pulsators, with av¬ 
erage ([Fe/H] 2 ) =—1.01 ±0.05. The three-parameter formula gives significantly 
lower value, ([Fe/HJa) =—1.21 ±0.07. This result is not unexpected. Since the 
three-parameter formula contains A 2 amplitude (Eq. 3), which is on average 
lower for BL stars {cf. Figs. 7 and 8 , Table 4), the calculated abundances are 
smaller. 

We have to deal with a question, which formula, two or three-parameter, 
gives proper results. It is a question, whether the lower values of A 2 ampli- 
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tudes for BL pulsators, are caused by their lower [Fe/H], or are the consequence 
of the Blazhko phenomenon, leading to the change of the fundamental mode 
shape. There are two arguments in favour of second explanation. First, the 
two-parameter [Fe/H] estimator (which has the higher dispersion only), doesn’t 
lead to the lower abundance of BL stars. Second, as it is visible on Fig. 8 , the 
interrelations between Fourier amplitudes of BL pulsators are disturbed. Result¬ 
ing conclusion is, that lower values of amplitude A 2 , are not directly connected 
with metallicity, but are rather a dynamic effect. In the case of BL pulsators we 
should use only the two-parameter estimator of [Fe/Hj. We have already shown, 
that using the two-parameter (P, (/ 31 ) estimator in case of the Flight curves, 
leads to metallicities consistent with spectroscopic values. We will return to the 
differences between our [Fe/H] estimators, while analyzing the LMC sample. 

No matter which formula is used, the Blazhko variables do not seem to prefer 
higher metallicities as suggested by Moskalik and Poretti (2003). 

4.2 The LMC pulsators 

Metallicity distributions of the LMC RRab stars are presented on Fig. 9. Un¬ 
expected result is visible for RRO-S stars metallicity distributions. In the case 
of the Galactic Bulge, both two and three-parameter [Fe/H] estimators lead 
to consistent results. In the LMC case, the distribution based on a three- 
parameter [Fe/H] estimator is broader and has a lower average value, ([Fe/H] 3 ) = 
-1.387±0.005, while ([Fe/Hja) =-1.218±0.004. The three-parameter [Fe/H] 
estimator leads to systematically lower metallicities, than the two-parameter 
one, with the difference increasing with decreasing [Fe/H]. 

Empirically, we may check which [Fe/H] estimator, two or three-parameter, 
gives proper results, by comparing with spectroscopic [Fe/H]. According to 
Clementini et al. (2003), average spectroscopic metallicity of 101 LMC RR Lyrae 
stars is —1.48 ±0.03 on the Harris (1996) metallicity scale. Individual [Fe/H] 
values for these stars are given in Gratton et al. (2004) who compared also 
spectroscopic [Fe/H] with photometric [Fe/H] for 19 RRab stars. Photometric 
metallicities were derived from V light curve and JK formula. It was found that 
the average difference between the JK photometric scale and the Harris scale is 
0.28±0.05. In the considered metallicity range ([Fe/H] > —1.8) JK photometric 
metallicities are on average higher by 0.05 ±0.03 than the HDS metallicities 
(c/. Section 2.1). Taking this two effects into account, we may state, that the 
HDS scale is more metal-rich than the Harris scale by 0.23 ±0.06. Converting 
metallicity given by Clementini et al. (2003), from Harris scale to HDS scale 
(used in this paper), we find —1.25±0.07. Comparing this value with average 
[Fe/H] derived from our formulae, we conclude that the value obtained with two- 
parameter formula is in very good agreement with spectroscopic [Fe/H], which is 
not the case for three-parameter formula. We conclude that the two-parameter 
formula is a better metallicity estimator for the LMC RRO-S stars. 

The difference between metallicities obtained from formulae with and with¬ 
out A 2 for the RRO-S LMC stars, not noticed in the case of Galactic Bulge, 
seems to be connected with physical properties of the stars of these systems. 
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Table 5: Averaged Fourier amplitudes and phases in three period bins, both for 
singly periodic and Blazhko variables from the LMC. The error of the mean is 
given in a parenthesis 



P<0.5 

Pe (0.5,0.6) 

P>0.6 


RRO-S 

RRO-BL 

RRO-S 

RRO-BL 

RRO-S 

RRO-BL 

Al 

0.237(2) 

0.227(3) 

0.196(1) 

0.190(2) 

0.153(1) 

0.147(4) 

A 2 

0.114(1) 

0.097(2) 

0.094(1) 

0.081(1) 

0.071(1) 

0.064(2) 

A 3 

0.081(1) 

0.056(2) 

0.0670(4) 

0.051(1) 

0.045(1) 

0.040(1) 

A 4 

0.056(1) 

0.037(1) 

0.0452(3) 

0.034(1) 

0.0314(4) 

0.028(1) 

A 5 

0.0410(4) 

0.031(2) 

0.0323(2) 

0.027(1) 

0.0232(4) 

0.021(1) 

021 

2.541(6) 

2.55(1) 

2.698(3) 

2.689(7) 

2.975(5) 

2.93(1) 

031 

5.27(1) 

5.18(3) 

5.613(7) 

5.50(1) 

6.15(1) 

6.05(3) 

041 

1.76(1) 

1.66(5) 

2.25(1) 

2.09(3) 

2.98(2) 

2.74(6) 

051 

4.50(2) 

4.50(6) 

5.13(1) 

4.99(4) 

5.82(3) 

5.61(9) 


The comparison of pulsational properties of stars, being members of our cali¬ 
brating sample, the Galactic Bulge and the LMC, is done, by construction of A 2 
vs. Ai plot - Fig. 10. For clarity, on the left panel we compare the calibrating 
and the Galactic Bulge samples, while on the right panel the Galactic Bulge 
and the LMC samples. Stars of the calibrating sample and the Galactic Bulge 
follow the same progression, while progressions of the LMC stars and the Galac¬ 
tic Bulge stars are inclined to each other. The effect is small, but visible, and 
is most pronounced for A 2 (Ai) relation. For this reason we should not use the 
three-parameter [Fe/H] estimator to the LMC stars. Only the formula without 
amplitude A 2 , will be used in further considerations. 

It appears that photometric method, is not at all universal, but slightly 
depends on the considered stellar system. This was not previously noticed, 
simply because the JK formula, doesn’t include amplitudes, and it seems, that 
they are not necessary, when one uses the V light curves. On the base of this 
study, we may state, that amplitudes are probably important in the case of 
/ light curves, however a more detailed calibration, based on a larger, more 
homogeneous sample is required. Special care should be taken, while applying 
the method to different stellar systems. 

Fourier amplitudes and phases vs. period, for singly periodic and Blazhko 
variables, are presented on Figs. 11 and 12. For clarity only one third of the 
stars are plotted. Similarly as for the Galactic Bulge stars, amplitudes of RRO- 
BL stars are smaller than for the RRO-S stars. Average Fourier parameters are 
presented in Table 5. The phases of both RRO-BL and RRO-S pulsators (Fig. 12), 
seem to follow the same sequence, however it is clear from Table 5, that indeed 
phases are also smaller for Blazhko pulsators. The smallest difference we note 
for the low order parameters, A 2 and (/) 2 i ■ Similar results were found by Alcock 
et al. (2003), for their MACHO LMC sample. 
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It is also worth noticing, that Blazhko pulsators prefer shorter periods. Only 
7 Blazhko stars have period longer than 0.7 days. As can be seen from period 
distribution, presented on Fig. 13, Blazhko variables with periods shorter than 
0.6 days, constitute 17% of all RRab stars. For longer periods, there is a signif¬ 
icant decrease of Blazhko effect incidence rate. Similar tendency is seen in the 
Galactic Bulge. 

Using the two-parameter [Fe/H] estimator, we can now compare the metallic- 
ity distributions of RRO-S and RRO-BL stars (left panel of Fig. 9). The average 
[Fe/H] values are —1.218 ±0.004 and —1.28 ±0.01 respectively, so they are sim¬ 
ilar, with average being lower for BL variables. The metallicity distribution 
for RRO-S stars is asymmetric, with a metal-rich tail, while the distribution for 
RRO-BL stars is approximately symmetric. 

Lower average metallicity for RRO-BL stars and a metal-rich tail for RRO-S 
stars, suggests, that Blazhko variables may prefer lower metallicities, oppo¬ 
site to Moskalik and Poretti (2003) suggestion. To verify this hypothesis we 
divide all the stars into three groups with different metallicities, and in each 
group we calculate the incidence rate for RRO-BL stars. Results are collected 
in Tables 6A-C. To avoid possible dependence of result, on the way groups are 
defined, division is made in three different manners: Table 6A - groups with 
approximately equal number of RRO-BL stars. Table 6B - groups with approxi¬ 
mately equal number of RRO-S stars and Table 6C - central group concentrated 
on the average metallicity of the RRO-S stars, with the width 0.4, and remaining 
two groups corresponding to the wings of the metallicity distributions. 

Regardless of the division scheme, the RRO-BL stars are most frequent 
(19 — 20%) in the lowest metallicity group, with [Fe/H] lower than approxi¬ 
mately — 1.4. In the remaining two groups of higher metallicity ([Fe/H] >« — 1.4) 
incidence rates are approximately the same (12 — 14%). 

Such an effect is not visible in the case of the Galactic Bulge stars. Here, due 
to the less numerous sample, we divide the stars into two groups, with metal¬ 
licities greater and lower than the average [Fe/H] for RRO-S stars. Obtained 
incidence rates of Blazhko variables are similar: 26 ±6% in the low metallicity 
group and 30±6% in high metallicity group. Taking into account the metallicity 
spread in the Galactic Bulge, ~ (—1.4, —0.8), this result is consistent with the 
LMG result, in the sense that for [Fe/H] higher than «—1.4, the incidence rate 
for BL stars is approximately constant. 

We have also checked the relations between [Fe/H] and dynamical parameters 
characterizing the Blazhko phenomenon, namely the period of the Blazhko cycle, 
Pb , and the amplitude of the highest additional frequency componenet, Ab- 
No clear dependence was found. 

4.3 Luminosity — metallicity relation for the LMC RRO-S 
pulsators 

Having determined metallicity for thousands of RR Lyrae stars we are able 
to derive the luminosity-metallicity relation with very high statistical accuracy. 
This relation is often assumed to be linear of the following form: My=a [Fe/H] ± 
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(3, however there is no agreement on the exact values of the a and (3 parameters. 
Derived values of the slope fall in a range 0.2 — 0.3 (see Olech et al. 2003 and 
references therein), while zero point estimates fall in a range 0.5 —0.7 at [Fe/H] = 
— 1.5 with the average value from ten different methods equal 0.59 ±0.03 (see 
Cacciari and Clementini 2003 for review). 

The Vo — [Fe/H] relation for singly periodic LMC RRab stars is presented 
on Fig. 14. Vo values are extinction free magnitudes directly taken from OGLE 
catalog (Soszyhski et al. 2003). We have omitted stars from OGLE fields SG17- 
SC20. They have on average lower magnitudes and form a more distant part 
of the warped LMG bar (Subramaniam 2003). For stars being overlaped on 
different fields, independent Vq measurments are averaged. Linear fit to the 
data (3397 RRab stars) leads to relation: 

Vo= 0.30[Fe/H] ±19.33 . . 

±0.02 ±0.02 

No evidence of nonlinearity is found. Obtained slope is rather steep, and agrees 
well with Olech et al. (2003) value (0.26±0.08), but is higher than values given 
by Glementini et al. (2003) (0.21 ±0.05) or Fernley et al. (1998) (0.18±0.03). 

To compare our estimation of zero point with those collected in Gacciari 
and Clementini (2003), we note, that their value of [Fe/H] = — 1.5 corresponds to 
AS'=7.03 (Fernley and Barnes 1997) which corresponds to [Fe/H]=-1.31 on HDS 
scale (Jurcsik 1995). Our relaton leads to Vo = 18.939±0.005 at [Fe/H] =—1.31. 
The absolute magnitude value depends on the assumed distance modulus to the 
LMC, /iLMC) but there is a large spread in the estimated values (see Clementini 
et al. 2003, specially their Fig. 8). Adopting the preferred value used by HST 
Key Project, /tlmc = 18.5±0.1, (Freedman et al. 2001) we obtain 0.44±0.10 
mag for the zero point in absolute magnitude (at [Fe/H] = — 1.31), brighter than 
the value given by Cacciari and Clementini (2003) (0.59±0.03). The difference 
may be due to overestimated OGLE reddenings given in Udalski et al. (1999) 
and adopted in Soszyhski et al. (2003) catalog, as noticed in the literature 
{e.g., Clementini et al. 2003). The average OGLE reddening is: E(i?—E)=0.143, 
to be compared with values adopted by HST Key Project (Freedman et al. 2001) 
and Fouque et al. (2003), E{B—V)=0.1, value given by Laney and Stobie (1994), 
E{B — V) =0.07, or Subramaniam (2005), E{B — V) = 0.06 (value transformed 
from E{V — I) using Fouque et al. 2003 coefficients). We note that decrease of 
the average OGLE reddening from E{B — V) = 0.143 to E{B — V) = 0.10 leads to 
decrease of absorption from Ay=0.47 to Ay=0.33 (assuming relation of Fouque 
et al. (2003): Ay = 3.3E(i? —K)) and hence shifts the zero point towards fainter 
magnitudes by 0.14. This leads to a very good agreement with Cacciari and 
Clementini (2003) value. 


5 Conclusions 

In order to examine the large data samples of the OGLE database: the Galactic 
Bulge and the LMC RR Lyrae samples, we have calibrated to the /-band the 
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Table 6 : The incidence rates of Blazhko variables in the LMC. The three Tables, 
A, B and C, correspond to the different manners of division of stars into groups, 
as described in text. 



[Fe/H] 

Nbl 

Ns 

% 

A: 


<-1.37 

235 

989 

19.2 ±1.1% 


(-1.37, -1.19) 

236 

1460 

13.9 ±0.8% 


>-1.19 

231 

1541 

13.0 ±0.8% 

B: 


<-1.32 

308 

1352 

18.6 ± 1 . 0 % 


(-1.32, -1.16) 

197 

1310 

13.1 ±0.9% 


>-1.16 

197 

1328 

12.9 ±0.9% 

C: 


<-1.42 

183 

702 

20.7 ±1.4% 


(-1.42, -1.02) 

436 

2644 

14.2 ±0.6% 


>- 1.02 

83 

644 

11.4±1.2% 


photometric method of determining the [Fe/H] of the RRab star, proposed by 
Kovacs and Zsoldos (1995). Our calibrating sample consist of 28 singly periodic, 
galactic RRab stars. The best two-parameter formula contains the period and 
the Fourier phase, ^ 31 , just as in case of using Flight curves. However the three- 
parameter formula, containing also the amplitude A 2 , was found to be a better 
[Fe/H] estimator, at least in the Galaxy. In case of the LMC, only the two- 
parameter formula gives the values of [Fe/H] which are consistent with direct 
spectroscopic determinations. It is a consequence of slightly different pulsational 
properties of the Galactic and the LMC stars. Only the two-parameter formula 
should be used for the LMC stars. 

Average metallicity of the Galactic Bulge RRO-S stars is —1.04±0.03, while 
for the LMC it is lower, —1.218±0.004. In both populations spread of individual 
metallicity values is larger than the error of the method, indicating intrinsic 
metallicity spread of these populations (of the order of 0.16 for the Galactic 
Bulge and 0.19 for the LMC). 

We showed, that the V photometric method, may be used to obtain reliable 
values of Blazhko variables metallicity, if the Fourier parameters of the light 
curve, prewithened with frequencies of additional componenets, are used. We 
assume that it also holds for the /-band, used thorough this paper. Also here, 
only the two-parameter formula should be used, to determine the metallicity of 
the Blazhko variables. It is a consequence of Blazhko phenomenon dynamics, 
that leads to the change of amplitudes of harmonics of the fundamental mode. 
Average amplitudes are smaller than for the singly periodic variables, both in 
the Galactic Bulge and in the LMC. 

In the Galactic Bulge, the incidence rate of Blazhko effect, doesn’t depend on 
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metallicity, and is equal ~25 —30% in the whole metallicity range, (—1.4,—0.8). 
In the LMC, Blazhko variables prefer lower metallicity values. For [Fe/H] < —1.4 
the incidence rate is significantly higher (~20%), while for higher metallicity 
values, it is approximately constant (~13%). Thus, the higher incidence rate 
of Blazhko variables in the Galactic Bulge, as compared to the LMC, cannot be 
attributed to the metallicity difference between the two stellar systems, as was 
suggested by Moskalik and Poretti (2003). 

Considering the global properties of Blazhko and singly periodic variables, 
we note that Blazhko pulsators prefer shorter periods, with strongly decreasing 
incidence rate for P>0.6 days. Fourier phases and harmonics of the amplitudes 
are lower, than for the singly periodic variables, as was also noticed by Alcock 
et al. (2003) for their MACHO LMC sample. The relations among Fourier 
amplitudes are disturbed in comparison with singly periodic variables. 

Estimation of metallicity for thousands of RR Lyrae stars allowed us to in¬ 
vestigate the luminosity-metallicity relation with very high statistical precision. 
We found a steep dependence of My on [Fe/H], with the slope equal 0.30±0.02. 
Our estimation of zero point of this reltion agrees with the literature values, if 
we assume that OGLE E{B — V) reddenings are overestimated by ~0.04 mag. 
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ing the Zangwill’s minimisation algorithm. 


6 APPENDIX: Analysis of the Blazhko variables 
from the LMC 

We describe each Eourier series representing the Blazhko variable light curve, 
with the following sequence: N = {No, A_, N+, Aa), hereafter shortened to N. 
Nq denotes a number of primary frequency harmonics increased by one, A_/_|_ 
stands for the number of harmonics accompanied by additional lower or higher 
secondary frequency, also increased by one and Aa is equal 1 if a frequency | A/| 
is detectable, or 0 otherwise. Only terms with A/(t^> 4 are considered. In order 
to find a sequence A that represents the Fourier sum, fitting the data best, the 
automatic procedure was developed. Conditions described below, were chosen, 
so that the algorithm works most effectively for the data for 30 pulsators of the 
Galactic Bulge sample, being strictly of BLl or BL2 type. This test is described 
at the end of the Appendix. The procedure can be described in the following 
steps: 

1 The OGLE frequencies of fundamental mode, fO, and of additional com¬ 
ponent of the highest amplitude, f 1, are used as an input data. Frequency 
separation, Df=f 1-fO, is calculated. 
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2 In the main loop of the procedure, all sensible N are fitted to the data. 
Following indices goes through: n0=3,7; np=0 ,nO+l; nm=0,n0+l; nD=0,1; 
ins=0,l. nO corresponds to Nq. If n0=7 is found, also higher values are 
checked, np is a number of components of the form i-fO—Df , while nm is 
a number of components of the form i-fO+Df . Depending on the Df sign 
they correspond to or N-. nD corresponds to N^. ins is equal 1 if 
the instrumental term of period 50,000 days is considered or 0 otherwise. 

3 All the bad data points rejected in the previous steps are recovered. For 
each N being tested, bad data points are rejected independently. 

4 The multiperiodic Fourier series with the frequencies determined by N, is 
fitted to the data. Dispersion of the fit, sig, is calculated. 

5 Bad data points (threshold 4sig) are rejected using following criterions: 

- First bad data point is always rejected, then procedure goes to step 4. 

- If at least one bad data point was already rejected, the next one is 
rejected only if np=0. If np^fO the point is kept in the sample. Please, 
notice that when np=0 and nniT^O we check for the frequency in the location 
it is expected (the Df sign gives us information about location of the 
additional frequency - on the right or left of the primary frequency). This 
is the way we protect the procedure, from finding BLI type with additional 
componenet on the wrong side of primary frequency, by rejecting to many 
data points. 

- The above condition doesn’t apply, if the correct Nq sequence was found 
(in the second loop of the algorithm). 

6 It is checked whether for all the terms, A/(t^> 4. If not, the procedure goes 
to step 2 and the next sequence N is considered. If all the terms satisfy 
the above condition, it is checked whether their number is greater by 2 
or more, from the number of frequencies of the previously remembered 
sequence. If not, the sequence is remembered, otherwise not. It protects 
the procedure from finding apparently good solutions by fitting the Fourier 
series of to high order. 

7 From all the remembered sequences the one with the smallest sig, Nq, 
is chosen. The procedure is started again, but only the frequencies de¬ 
termined by the Nq are fitted. There is no limit on the bad data point 
rejection and the initial frequencies are adjusted using the Zangwill (1969) 
conjugate directions procedure. 

8 Procedure goes to step 1 and the next star is examined. 

The algorithm was tested (without frequency adjustment) on the 30 Galactic 
Bulge, strictly BLI or BL2 pulsators, described by Moskalik and Poretti (2003) 
by simple N. As an input we used the primary frequency and one secondary 
frequency of the highest amplitude. Only in two cases different sequences were 
found. For BW10_V44 the (4, 0, 3, 1) sequence was found, instead of (4, 0, 
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4, 0), but it didn’t lead to significant change of Fourier parameters. The only 
serious problem was accounted for BW10_V41, were the BL2 type (4, 3, 3, 0) 
was found, instead of a BLl, (4, 0, 3, 0), which led to the significant change in 
the Fourier parameters. 

Based on these tests we conclude, that in vast majority of cases, our algo¬ 
rithm is satisfactory to obtain reliable Fourier parameters for the fundamental 
mode. 
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Figure 1: Spectroscopic vs. photometric [Fe/H] estimated with two (left panel) 
and three-parameter calibration (right). 


Am -0.6 

■ a) 

-0.6 

■ b) 

-0.4 

: ..... 

-0.4 

K K 


i 'y t *»••• 


: > : > 

-0.2 

;■ 1. ;■ i . A - 

-0.2 

t i 

■ \ ■ 'V 





0 


0 


0.2 


0.2 

S ; 

- • 

0.4 

... 1 ... 1 ... 1 ... 1 .. . 

0.4 

... 1 ... 1 ... 1 ... 1 .. . 


0 0.4 0.8 1.2 1.6 2 0 0.4 0.8 1.2 1.6 2 


Figure 2: The light curve of the RR0-BL2 pulsator, BW6_V7 (OGLE-1 sample); 
a) light curve phased with the fundamental mode frequency, b) light curve after 
the secondary frequencies have been filtered out. 






20 




[Fe/H]2 [Fe/H]3 

Figure 3: The [Fe/H] distributions obtained with two-parameter (left) and three- 
parameter formula (right) for RRO-S (top) and RRO-BL (bottom) pulsators. For 
each distribution the average abundance, ([Fe/H]), and the sample dispersion, 
a, is given. 
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Figure 4: The peak-to-peak amplitude. Am, vs. period for the Galactic Bulge 
RRO-S stars of different [Fe/H], estimated with two (left panel) or three (right 
panel) parameter formula. Stars with different [Fe/H] are plotted with different 
symbols: (□) [Fe/H] >-0.8; (x) -1.2 < [Fe/H] <-0.8; (v) [Fe/H] <-1.2. 
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Figure 5: The Fourier phases, (pzi and (/) 4 i, vs. period for the Galactic Bulge 
RRO-S stars. Abundances are estimated with two (left) or three-parameter 
formula (right). Line of [Fe/H] = —1 is plotted. Symbols are the same as in 
Fig. 4 
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Figure 6: The Fourier amplitude, Ai, vs. period for RRO-S (•) and RRO-BL (0) 
pulsators of the Galactic Bulge. 
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Figure 8: The Fourier amplitudes, A 2 —A 5 , vs. the Fourier amplitude of the 
fundamental mode, Ai, for the Galactic Bulge pulsators. Symbols are the same 
as in Fig. 6 
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Figure 9: The [Fe/H] distributions obtained with two-parameter (left) and three- 
parameter formula (right) for RRO-S (top) and RRO-BL (bottom) LMC pul¬ 
sators. For each distribution the average abundance, ([Fe/H]), and the sample 
dispersion, a, is given. 
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Figure 10: The comparison of A 2 vs. Ai relations for the calibrating sample, 
the Galactic Bulge and the LMC. For clarity one third of the LMC stars are 
plotted. 
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Figure 11: The Fourier amplitude, Ai vs. period, for RRO-S (•) and RRO-BL 
(0) pulsators of the LMC. 
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Figure 13: Distribution of periods for RRO-S and RRO-BL pulsators of the LMC. 
The Blazhko variables incidence rate is given for four period bins. 
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Figure 14: The extinction free apparent magnitude vs. metallicity relation for 
the LMC RRO-S pulsators. 
















